General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of som 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 


i 


i 


•-*r 


aiiiMiriAtfl 




- \ 




' 


. 

. *’ , 












. 








■ t . 


(NASA-CB- 1755<e2) THE IBFLIENCE OF TH3 
HELI CSPHEBIC CUBBEBT SHEET ANE ANGULAR 
SEPARATION ON FLARE ACCEL EBA 1 ED SCLAB WIKD 
(Stanford Oniv. ) 17 p HC A02/HE AO 1 

CSCL 03B G3/92 


* 7 * 


NB5-2S9C7 


Urc. la e 
24057 


f 





The Influence of the Heliospheric Current Sheet 
and Angular Separation oh Flare Accelerated Solar Wind 


Harald M. Henning, Philip H. Scbcrrer and J, Todd Hocksema 
Center for Space Science and Astrophysics 
Stanford University, Stanford, CA 


CSSA-ASTRO-85-I7 


April 1085 


Submitted to: Journal of Geophysical Research 


Office of Naval Research Contract NQ0014-76-C-0207 


National Aeronautics and Space Administration Grant 



GR5-P20-559 


National Science Foundation Grant ATM77-20580 




O 


© 


c 


D 


O 


c 


o 


G 


Abstract. A complete set of major flares has beco used to Investigate the effect of the 

heliospheric current sheet on the magnitude of the flare-associated disturbance measured at 

Earth. It was found that disturbances associated with flares located on the same side of the 

current sheet as Earth were of larger magnitude than those associated with flares located 

such that the flare-acccicratcd material would have to cross the current sheet. It was also 

found that the angular separation between the flare position and Earth has a strong effect 

* 

on the magnitude of tin disturbance. A larger angular separation tended to result in a 
smaller disturbance. Thirdly, it was determined that flares tend to occur near the helios- 
pheric current sheet. 

1. Introduction 

» 

It is known that flares produce disturbances in the solar wind resulting in geomagnetic 
activity at Earth. However, the magnitude of such disturbances cannot yet be predicted for 
any given flare. This paper presents the results of a statistical analysis which separates and 
identifies some of the quantities which determine the relative magnitude of a flare-associated 
disturbance. The main thrust of the analysis is to determine the effect of the relationship 
between the location of the flare and the location of the heliospheric current sheet (HCS). To 
avoid bias in the estimate of an HCS effect, however, one must also examine any effect due to 
the position of the flare on the solar disk. This analysis produces two major results: First, that 
flares occurring on the opposite side of the HCS from tbt earth produce smaller solar wind and 
geomagnetic disturbances than those on the same side. Second, we find that the disturbance 
decreases with increasing angular distance of the flare from the 6ub-earlh point. 

Three kinds of data were used in the present study. The first is the list of major flares 
fcs defined by Dodson and Hcdcman (1071, 1075, 1081). The second is the location of the HCS 
as computed by Hockscma (1082, 1083, 1084, 1085). The third is the indicator of solar wind 
response to the flare, and has been taken to be the geomagnetic disturbance index Dst and 
aolar wind velocity when available. 
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I. Data 

The flares used in this analysis are taken from the UAG reports of the Comprehensive 
Flare fndex (CFI) compiled by Helen W. Dodson and E, Ruth Hedeman (1071, 1075, 1001), 
These records contain a complete list of all observed rm.jor flares from 1055-1070. The CFI is 
comprised of five components: 

1) Importance of ionizing radiation as indicated by time-associated Short Wave Fade 

* 

or Sudden Ionospheric Disturbance (Scale 0-3), 

2) Importance of Ha flare (Scale 0-3), 

3) Magnitude of f^lOcm flux; (characteristic of log of flux in units of 
lO~ n Wm~ 7 He~ l ); (Scale 0-4). 

4) Dynamic spectrum; (Type II «= 1, Continuum = 2, Type IV with duration > 10 
min = 3); (Scale 0-3). 

5) Magnitude of ca. 200Aj7/z flux; (characteristic of log of flux in units of 
10'” Wm' 2 (Scale 0-5). 

The total CFI of a flare is then Just the sum of the 5 components, with a possible range 
of 1 to 18. 

A major flare is defined as one which was well above average in either ionizing, II a, or 
radio frequency radiation (corresponding CFI component >3), or exhibited a Type II burst or 
Type IV burst with longer tlan 10 min duration (Dodson and Hedeman 1971, 1975, 1981), 
From this still formidable list the data for this analysis was selected by including only flares 
with a total CFI of 7 or greater. For the interval from 1968 through 1979 the data includes 
505 flares. Missing data for the flare position or solar wind parameters reduces this number to 
slightly more than 400. All these flares were used for the part of the analysis determining the 
disk-position dependence of solar wind disturbances. 

The location of the heliospheric current sheet (HCS) is derived from the potential field 
model calculations done at Stanford (Hocksema 1982, 1983, 1984, 1985). The HCS data starts 
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in the middle of 1070, and continues to the present. For the HCS analysis then,' the flares 
from 1070 to 1070 were used. The majority of these occurred in 1078 and 1070, as the solar 
cycle maximum was approaching, There are 181 flares in this interval, of which slightly more 
than HO remain after elimination of those with missing positions or solar wind data. 

The geomagnetic index D S ti which is continuously available and proportional to solar 
wind speed and D, (Murayama 1082), and the solar wind speed hourly averages observed by 
IS12E3 are used to measure the flare generated disturbance. Figures 1 and 2 show superposed 
epoch analyses of 6olar wind 6pccd and D sr about flare times. These plots show a strong peak 
in the disturbance of the solar wind and D$r 2 to 4 days after the flare. Note that a higher 
value for solar wind speed is the signature of a disturbance, whereas for Dst it is a lower, i.e. 
greater negative, value. We chose to examine 24-hour averages of these two indices to study 
the effect of flare position on the disk and relative to the HCS. Use of the 24-hour averages 
has the advantage that it avoids the subjectivity that occurs when one attempts to identify 
the source of a particular disturbance. An analysis of 24-hour intervals in the period from 1.0 
to 4.6 days after the flare showed that averages of the interval from 2.0 to 3.0 days after the 
flare demonstrate the effect of the HCS and flare position moat strongly. We found that there 
was a strong HCS effect evident when examining any 24-hour interval between 1.5 and 4.6 
days after the flare, beginning sharply at 1.5 days and a gradual decrease commencing 3 days 
after the flare and continuing past 4.5 days. Thus, for each flare in our study, a 24 hour aver- 
age was calculated for the interval starting 2 days after the flare. 

The time interval chosen allows for a propagation velocity of the flare-accelerated 
material ranging from 577 km/scc to 8G6 km/sec. Note that the 24 hour average of the solar 
wind speed itself will usually be substantially less than this, since the average also includes 
quiet solar wind. Very likely some of the flare associated disturbances may be missed, but, 
without making the time intervals too large and further diluting the site of a disturbance, this 
problem cannot be avoided in a completely objective analysis where we want to avoid the 
uncertainty in associating a specific solar wind or Dsx disturbance with a particular flare. 
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8. Analysis 
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8.1* The Heliospheric Current Sheet 
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To determine whether a Dare-associated disturbance had to cross the HCS to reach 

Earth we looked at the source surface calculated at Stanford with the potential field model 

(Hoekscma 1082, 1083, 1084, 1083). We assume the steady, undisturbed solar wind carries 

the magnetic structure of the source surface into interplanetary space, reaching Earth after 

an average of 5 days. The flare material is assumed to take an average of 2.5 days to 

reach Earth, Thus wc can estimate the magnetic structure at the flare position on the sun 

and at Earth when the disturbance arrives. Specifically, if the polarities of the magnetic 

Sold (cither inward directed or outward directed) are different, the flare is considered to be 

* 

on the opposite side of the HCS, l.e. the flare material has to cross the HCS to reach Earth. 
If the polarities are the same, the flare is considered to be on the Bame Bide of the HCS. If 
either of the positions was very close to the HCS, i.e. the Held magnitude was very small 
(less than 0.1 (iT on the source. surface), the flare was eliminated from the analysis. There 
were four such instances. In addition, two cases where the flare material would have pro- 
pagated tangential to the HCS were eliminated from the 6tudy, since the latitudinal extent 
of the HCS is uncertain by a few degrees. There were three cases where flare material had 
to cross two HCS’s. The final flare 6ct for the solar wind speed analysis contains 135 flares, 
of which 75 are on the same side of the HCS and 60 on the opposite side. There were 2 
more flares in each category available for the D ST analysis. 

We define a quantity D , which we call the difference measure, that provides an indi- 
cation of the significance of the effect the HCS has on flare- associated disturbances as fol- 
lows : 


far.... +3x „, ) ’ 

where X is the quantity with which we are measuring the disturbance, e.g. solar wind 
speed 3-4 days after the flare. 
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la words, D is the difference between the average quantity (solar wind speed, D S r) 
for flares on the same side off the HCS and the overage quantity for flares on the opposite 
side of the HCS, as measured in units of the average error of the mean. Thus, a large 
absolute value of D is an Indicator of an effect. 

Figure 3 shows the values of D for solar wind speed 2 to 3 days after the flare. 
Large positive D means a significantly larger value of average solar wind speed for flares 
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on the same side of the HCS than for flares on the opposite side, In panels a-e the flares 
were grouped depending on their value for different components of t**e CFI. The flares in 
figure 3a were grouped by the value of the first component of the CFI, i.c. by the impor- 
tance of ionizing radiation. The value of D for all flares which had ionizing radiation 
importance of 1 or greater is shown at the abscissa "1+”, the value of D for all flares with 
ionizing radiation importance of 2 or greater at abscissa ”2+”, and so on. Similarly, figure 
3b was grouped by Ha importance, 3c by the magnitude of 10cm flux, 3d by the dynamic 
spectrum, and 3c by the magnitude of 200MHz flux. Note that the first point in each 
panel includes all flares with a CFI component value of 0 or greater, i.c. all flares in the 
data. The reason for thi.i method of grouping is that 6incc flares arc selected on the basis 
of a value of total CFI >7, a value of 0 in one component of the CFI implies on the aver- 
age a high value in another component of the CFI, This bias prevents us from examining 
flares grouped by a particular value of a particular CFI component, rather we examine 
Bares grouped according to whether a CFI component exceeds a given threshold. The 
adopted grouping enables one to sec better the effect large flares have as opposed to smaller 
ones. 

In the final panel, f, flares were categorized by their total CFI value, in the method 
described above. 

Figure 4 has the same format as figure 3, but shows the response of Dst rather than 
the solar wind speed. In this figure, a large negative value of D indicates an effect, since a 
more negative value of the D st index indicates greater geomagnetic activity. 
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The first and most important thing to be learned from these plots is lhat'there is a 
definite HCS effect. The flares on the same side of the IICS as the earth result in a larger 
disturbance, i.e. greater solar wind speed and more negative D ST , The size of this effect is 
quite large, on the order of 4a (solar wind speed) or 5a ( D sr ), which corresponds to about 
i G5 km/scc difference in solar wind speed or 18 in D S r units. In order to check for con- 
sistency of the effect, the flares were divide'.! into two time intervals and the calculations 
done separately for each. The resulting plots were practically identical to each ‘other and 
to the calculations from the full data (figures 3 and 4), demonstrating the consistency con- 
vincingly. 

None of the components of CFI, nor the total CFI, seems to produce a significantly 
larger or smaller effect. All the plots look fairly similar, even though for each panel 
different flares determined the points for CFI component values greater than 0. Thus none 
of the components of CFI stands out as a determinant for whether a flare will cause a large 
disturbance or not. Nor docs total CFI provide a better organization, although it certainly 
provides a more complete list of major flares than any of its separate components. 

One is tempted to note that.in almost all panels, for a higher value of the CFI com- 
ponent the absolute value of D is getting smaller, i.e. perhaps for the truly strongest flares 
the HCS effect is less. However, another factor contributes to produce this decrease in D , 
This factor is the definition of D . For the points at larger values of the CFI component 
there are fewer flares, and thus the error of the mean s larger, giving a smaller value of D . 
If one examines the actual solar wind speed or D ST difference between the average for 
flares on the same side and the average for flares on the opposite side, the tendency is actu- 
ally for the difference to be a larger number. Therefore, contrary to the immediate 
assumption, it appears that the HCS effect is proportionally stronger for stronger flares. 
However, the uncertainty of the numbers, as demonstrated by the decrease of D in the 
plots, means that we cannot provide a concrete answer to this question about large flares. 
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Summarizing this flection, by use of these graphs we have successfully demonstrated 
that there exists on HCS effect, given an idea of its significance and determined that no 
component of CFI, i.e. no energy range, stands out as a strong organizer of solar wind dis- 
turbances. 

3.2. Angular Separation 

In the preceding analysis, we found that the position of the flare with respect to the 
HCS, i.e. same or opposite side, resulted in a significant difference in the size of the flare- 
associated disturbance. However, to get an unbiased estimate of an effect the HCS has on 
such disturbances one has to consider a possible systematic effect associated with the angu- 
lar distance of the flare from tbe sub-earth position, it is clear that the farther a flare is 
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separated in angular distance from the sub-earth point, the more likely it is that an HCS 
lies in-botween. Thus there is a likely systematic bias that disturbances which have to 
cross the HCS arc due to flares which arc at larger angular distances from Earth. We have 
examined the solar wind response at Earth to a flare as a function of the angular distance 
between the flare and the sub-earth point. In addition to later enabling us, in our main 
analysis, to separate the HCS effect from this distance effect, this investigation is, of 
course, interesting in its own right. 

As a first test of the systematic bias, wc found that indeed, for tbe period starting in 
1076, where the HCS position is known, and ending in 1079, the end of the CFI compila- 
tion, the average angular distance from Earth for flares on the Bame side of the HCS was 
42 * , while for those on the opposite side it was 67* . 

To judge the effect of this bias we examined tbe dependence of tbe solar wind 
response, as measured by magnitude of solar wind speed or D$t, on the angular separation 
between the flare and Earth. Since we are not restricted by having to know the location of 
the HCS, we have taken all flares in our data with known positions, starting with the year 
1063 and running through 1070. The angular distance to the sub>earth point at the aver- 
age time when flare material is assumed to reach Earth, i.e. after 2.6 days, was calculated 
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for each Dare, The flares were then binned by distance and the average solar wind speed 
and Dst plotted in figures 6a and fla respectively. Figures 6b and Ob show the results o* a 
similar calculation, except that only longitudinal separation was considered. 

The graphs clearly show nn e#ect on tbi average size of the disturbance. The farther 
Earth is separated in angular distance from the position of the flare, the weaker the distur- 
bance measured at Earth. It is debatable whether the dependence is a monotonic decrease 
of the disturbance with angular distance, or whether there is a region of about £40* 
around the flare position where the disturbance is large, and outside of this ’’bubble” the 
effect of the flare is much smaller (especially evident for Dst)- However, the scatter plot of 
all flares seems to support a monotonic decrease. 

Also, one can argue that there seems to be an East-West asymmetry. Fcrr positive 
longitudinal separation th >i average magnitude of a disturbance, as measured by high solar 
wind speed or large negative Dsn seems larger than for corresponding negative longitudi- 
nal separation. Thus, it suggests that a flare on the western hemisphere of the sun may 
have a larger effect than one on the eastern hemisphere. This is consistent with the report 
that streams from flares on the western hemisphere have a speed 60% higher than those 
from eastern flares (Pudovkin et al 1079). However, since the uncertainty is fairly large we 
will consider only the angular distance in the further analysis. 

If we assume a roughly linear monotonic decrease, then wc can estimate the size of 
the systematic effect introduced into our HCS analysis. Wc calculated the average angular 
distance for flares on the same side of the HCS, 42* , and the opposite side, 67* . Figure 5a 
shows a difference of about 15 km/6cc in average solar wind speed between these two 
values. Figure 6a provides the same information for Dsr > a difference of about 3 D$t 
units. 

8.3. The HCS Effect 

Returning to our main investigation of the effect of the HCS on the disturbances 
caused by flares, wc must try to remove the effect of angular distance. To do so, a straight 
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tine was fitted to & plot of disturbance (as measured by solar wind speed or D ST ' measured 
at Garth 2.0 to 3.0 days after the flare) vs. angular separation for alt flares used in the 
HCS analysis (figures 7a and 7b). These lines then represents the effect of angular dis- 
tance! although one cannot totally separate the effects since at larger angular distance 
there are more flares on the opposite side of the HCS, and vice versa for smaller angular 
distance, influencing the slope of the line such that the HCS effect is weakened by an unk- 
nown, although certainly little, amount. In the next step the flares were separated accord- 
ing to whether they were on the same or opposite side of the HCS. For each flare the dis- 
tance from the line (in appropriate units, i.e. km/scc for solar wind speed, or D ST units) 
was calculated. Each group, i.e. same and opposite side flares, has a roughly Gaussian dis- 
tribution. The centroid and error of the mean of these distributions were calculated for 

4 

both groups of flares. Then thi difference between the centroid for flares on the same side 
of the HCS and the ccntioid for flares on the opposite side was taken, producing the fol- 
lowing results : 

i* 

Solar Wind Speed : some -opf. <= , 13.4±10.5 km / tee (98.7%>0) 

Dst * fame -opp. =-15.‘i±‘1.0 (00.0% <0) 

Referring to the values quoted in section 3.1 and 3.2, one secs that the angular dis- 
tance effect accounts for roughly one fourth of the difference demonstrated in figures 1 and 
2. The other three fourths result from the HCS effect. 

3.4, Distance from the HCS 

As an aside resulting from this analysis, the distance of each flare from the nearest 
HCS wa3 calculated, giving an average of 16* for flares from September of 1077 to the end 
of 1070, a period of continuous high activity. There was no significant difference between 
the average distance for flares on the same side of the HCS and the average distance for 
those on the opposite side. A set of artificial flare positions was then determined by keep- 
ing the actual latitudes (to insure a realistic distribution in latitude) but assigning random 
Carrington longitudes within this interval. Then the average distance to the HCS for this 
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act was calculated. This process was repeated 400 times. In none of these 400 trials did 
the average come as low as the true average, the lowest being only 17.5* . The overall 
average distance between random positions and the IICS for all iterations was 21*. This 
average is significantly greater than the average resulting from the true flare position. This 
confirms the tendency for flares to occur near the IICS (Dittmer 1076). 

4» Conclusion , 

In this analysis we have investigated the effect of angular separation and the HCS on 
flare associated disturbances measured at Earth. We have found that there is a clear associa- 
tion between proximity of the Dare to the sub-earth point and stronger disturbances. This 
relation accounts for a fourth of the difference noted between disturbances from flares on the 
same and opposite sides of the HCS. A strong effect associated solely with the HCS remains at 
a greater than 08% confidence level. Flares on the same side of the HCS tend to produce 
larger disturbances than flares on the opposite side. This influence of the HCS may be 
explained by the fact that solar wind speed tends to be at a minimum along the HCS (Dorrini 
1081), and thus disturbances in propagating to Earth woutd have to interact with this slower 
plasma, perhaps weakening them. 
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Flg»3. Graphs of the difference measure 
for solar wind speed vs. magnitude of t 
flare as measured with a) importance 
Ionizing radiation, b) importance of 7 
flare, c) magnitude of 10cm flux, 
dynamic spectrum, e) magnitude of 2007: 
flux, and f) total CFI. in each panel, 
runs from 0 (bottom) to 7 (top). 
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Flg.4. The same as figure 3, but for D ST 
rather than solar wind speed. The magni- 
tude of D is greater than for solar wind 
speed, otherwise the graphs are very simi- 
lar. In each panel, D runs from -7 (bot- 
tom) to 0 (top). 
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Flg.7. Graphs of a) solar wind speed and b) D sr 2 to 3 days after a flare vs. angular dis- 
tance. The line is a least squares fit to all points, and represents the effect of angular dis- 
tance on the magnitude of the flare disturbance. Flares occuring on the opposite side of the 
HCS from Earth arc marked with an o , those on the same 6idc with an x. 




